Toward a comprehensive language for biological systems

Rule-based modeling has become a powerful approach for modeling intracellular networks, which are characterized by rich molecular diversity. Truly comprehensive models of cell behavior, however, must address spatial complexity at both the intracellular level and at the level of interacting populations of cells, and will require richer modeling languages and tools. A recent paper in BMC Systems Biology represents a signifcant step toward the development of a unified modeling language and software platform for the development of multi-level, multiscale biological models

Computer Simulation of Cellular Patterning Within the Drosophila Pupal Eye

We present a computer simulation and associated experimental validation of assembly of glial-like support cells into the interweaving hexagonal lattice that spans the Drosophila pupal eye. This process of cell movements organizes the ommatidial array into a functional pattern. Unlike earlier simulations that focused on the arrangements of cells within individual ommatidia, here we examine the local movements that lead to large-scale organization of the emerging eye field. Simulations based on our experimental observations of cell adhesion, cell death, and cell movement successfully patterned a tracing of an emerging wild-type pupal eye. Surprisingly, altering cell adhesion had only a mild effect on patterning, contradicting our previous hypothesis that the patterning was primarily the result of preferential adhesion between IRM-class surface proteins. Instead, our simulations highlighted the importance of programmed cell death (PCD) as well as a previously unappreciated variable: the expansion of cells' apical surface areas, which promoted rearrangement of neighboring cells. We tested this prediction experimentally by preventing expansion in the apical area of individual cells: patterning was disrupted in a manner predicted by our simulations. Our work demonstrates the value of combining computer simulation with in vivo experiments to uncover novel mechanisms that are perpetuated throughout the eye field. It also demonstrates the utility of the Glazier–Graner–Hogeweg model (GGH) for modeling the links between local cellular interactions and emergent properties of developing epithelia as well as predicting unanticipated results in vivo

MDCK Cystogenesis Driven by Cell Stabilization within Computational Analogues

The study of epithelial morphogenesis is fundamental to increasing our understanding of organ function and disease. Great progress has been made through study of culture systems such as Madin-Darby canine kidney (MDCK) cells, but many aspects of even simple morphogenesis remain unclear. For example, are specific cell actions tightly coupled to the characteristics of the cell's environment or are they more often cell state dependent? How does the single lumen, single cell layer cyst consistently emerge from a variety of cell actions? To improve insight, we instantiated in silico analogues that used hypothesized cell behavior mechanisms to mimic MDCK cystogenesis. We tested them through in vitro experimentation and quantitative validation. We observed novel growth patterns, including a cell behavior shift that began around day five of growth. We created agent-oriented analogues that used the cellular Potts model along with an Iterative Refinement protocol. Following several refinements, we achieved a degree of validation for two separate mechanisms. Both survived falsification and achieved prespecified measures of similarity to cell culture properties. In silico components and mechanisms mapped to in vitro counterparts. In silico, the axis of cell division significantly affects lumen number without changing cell number or cyst size. Reducing the amount of in silico luminal cell death had limited effect on cystogenesis. Simulations provide an observable theory for cystogenesis based on hypothesized, cell-level operating principles

3D Multi-Cell Simulation of Tumor Growth and Angiogenesis

We present a 3D multi-cell simulation of a generic simplification of vascular tumor growth which can be easily extended and adapted to describe more specific vascular tumor types and host tissues. Initially, tumor cells proliferate as they take up the oxygen which the pre-existing vasculature supplies. The tumor grows exponentially. When the oxygen level drops below a threshold, the tumor cells become hypoxic and start secreting pro-angiogenic factors. At this stage, the tumor reaches a maximum diameter characteristic of an avascular tumor spheroid. The endothelial cells in the pre-existing vasculature respond to the pro-angiogenic factors both by chemotaxing towards higher concentrations of pro-angiogenic factors and by forming new blood vessels via angiogenesis. The tumor-induced vasculature increases the growth rate of the resulting vascularized solid tumor compared to an avascular tumor, allowing the tumor to grow beyond the spheroid in these linear-growth phases. First, in the linear-spherical phase of growth, the tumor remains spherical while its volume increases. Second, in the linear-cylindrical phase of growth the tumor elongates into a cylinder. Finally, in the linear-sheet phase of growth, tumor growth accelerates as the tumor changes from cylindrical to paddle-shaped. Substantial periods during which the tumor grows slowly or not at all separate the exponential from the linear-spherical and the linear-spherical from the linear-cylindrical growth phases. In contrast to other simulations in which avascular tumors remain spherical, our simulated avascular tumors form cylinders following the blood vessels, leading to a different distribution of hypoxic cells within the tumor. Our simulations cover time periods which are long enough to produce a range of biologically reasonable complex morphologies, allowing us to study how tumor-induced angiogenesis affects the growth rate, size and morphology of simulated tumors

Adhesion Failures Determine the Pattern of Choroidal Neovascularization in the Eye: A Computer Simulation Study

Choroidal neovascularization (CNV) of the macular area of the retina is the major cause of severe vision loss in adults. In CNV, after choriocapillaries initially penetrate Bruch's membrane (BrM), invading vessels may regress or expand (CNV initiation). Next, during Early and Late CNV, the expanding vasculature usually spreads in one of three distinct patterns: in a layer between BrM and the retinal pigment epithelium (sub-RPE or Type 1 CNV), in a layer between the RPE and the photoreceptors (sub-retinal or Type 2 CNV) or in both loci simultaneously (combined pattern or Type 3 CNV). While most studies hypothesize that CNV primarily results from growth-factor effects or holes in BrM, our three-dimensional simulations of multi-cell model of the normal and pathological maculae recapitulate the three growth patterns, under the hypothesis that CNV results from combinations of impairment of: 1) RPE-RPE epithelial junctional adhesion, 2) Adhesion of the RPE basement membrane complex to BrM (RPE-BrM adhesion), and 3) Adhesion of the RPE to the photoreceptor outer segments (RPE-POS adhesion). Our key findings are that when an endothelial tip cell penetrates BrM: 1) RPE with normal epithelial junctions, basal attachment to BrM and apical attachment to POS resists CNV. 2) Small holes in BrM do not, by themselves, initiate CNV. 3) RPE with normal epithelial junctions and normal apical RPE-POS adhesion, but weak adhesion to BrM (e.g. due to lipid accumulation in BrM) results in Early sub-RPE CNV. 4) Normal adhesion of RBaM to BrM, but reduced apical RPE-POS or epithelial RPE-RPE adhesion (e.g. due to inflammation) results in Early sub-retinal CNV. 5) Simultaneous reduction in RPE-RPE epithelial binding and RPE-BrM adhesion results in either sub-RPE or sub-retinal CNV which often progresses to combined pattern CNV. These findings suggest that defects in adhesion dominate CNV initiation and progression

Comparing individual-based approaches to modelling the self-organization of multicellular tissues.

The coordinated behaviour of populations of cells plays a central role in tissue growth and renewal. Cells react to their microenvironment by modulating processes such as movement, growth and proliferation, and signalling. Alongside experimental studies, computational models offer a useful means by which to investigate these processes. To this end a variety of cell-based modelling approaches have been developed, ranging from lattice-based cellular automata to lattice-free models that treat cells as point-like particles or extended shapes. However, it remains unclear how these approaches compare when applied to the same biological problem, and what differences in behaviour are due to different model assumptions and abstractions. Here, we exploit the availability of an implementation of five popular cell-based modelling approaches within a consistent computational framework, Chaste (http://www.cs.ox.ac.uk/chaste). This framework allows one to easily change constitutive assumptions within these models. In each case we provide full details of all technical aspects of our model implementations. We compare model implementations using four case studies, chosen to reflect the key cellular processes of proliferation, adhesion, and short- and long-range signalling. These case studies demonstrate the applicability of each model and provide a guide for model usage

A Multi-cell, Multi-scale Model of Vertebrate Segmentation and Somite Formation

Somitogenesis, the formation of the body's primary segmental structure common to all vertebrate development, requires coordination between biological mechanisms at several scales. Explaining how these mechanisms interact across scales and how events are coordinated in space and time is necessary for a complete understanding of somitogenesis and its evolutionary flexibility. So far, mechanisms of somitogenesis have been studied independently. To test the consistency, integrability and combined explanatory power of current prevailing hypotheses, we built an integrated clock-and-wavefront model including submodels of the intracellular segmentation clock, intercellular segmentation-clock coupling via Delta/Notch signaling, an FGF8 determination front, delayed differentiation, clock-wavefront readout, and differential-cell-cell-adhesion-driven cell sorting. We identify inconsistencies between existing submodels and gaps in the current understanding of somitogenesis mechanisms, and propose novel submodels and extensions of existing submodels where necessary. For reasonable initial conditions, 2D simulations of our model robustly generate spatially and temporally regular somites, realistic dynamic morphologies and spontaneous emergence of anterior-traveling stripes of Lfng. We show that these traveling stripes are pseudo-waves rather than true propagating waves. Our model is flexible enough to generate interspecies-like variation in somite size in response to changes in the PSM growth rate and segmentation-clock period, and in the number and width of Lfng stripes in response to changes in the PSM growth rate, segmentation-clock period and PSM length

A cell-based computational model of early embryogenesis coupling mechanical behaviour and gene regulation

The study of multicellular development is grounded in two complementary domains: cell biomechanics, which examines how physical forces shape the embryo, and genetic regulation and molecular signalling, which concern how cells determine their states and behaviours. Integrating both sides into a unified framework is crucial to fully understand the self-organized dynamics of morphogenesis. Here we introduce MecaGen, an integrative modelling platform enabling the hypothesis-driven simulation of these dual processes via the coupling between mechanical and chemical variables. Our approach relies upon a minimal ‘cell behaviour ontology’ comprising mesenchymal and epithelial cells and their associated behaviours. MecaGen enables the specification and control of complex collective movements in 3D space through a biologically relevant gene regulatory network and parameter space exploration. Three case studies investigating pattern formation, epithelial differentiation and tissue tectonics in zebrafish early embryogenesis, the latter with quantitative comparison to live imaging data, demonstrate the validity and usefulness of our framework

A Neutrophil Phenotype Model for Extracorporeal Treatment of Sepsis

Neutrophils play a central role in eliminating bacterial pathogens, but may also contribute to end-organ damage in sepsis. Interleukin-8 (IL-8), a key modulator of neutrophil function, signals through neutrophil specific surface receptors CXCR-1 and CXCR-2. In this study a mechanistic computational model was used to evaluate and deploy an extracorporeal sepsis treatment which modulates CXCR-1/2 levels. First, a simplified mechanistic computational model of IL-8 mediated activation of CXCR-1/2 receptors was developed, containing 16 ODEs and 43 parameters. Receptor level dynamics and systemic parameters were coupled with multiple neutrophil phenotypes to generate dynamic populations of activated neutrophils which reduce pathogen load, and/or primed neutrophils which cause adverse tissue damage when misdirected. The mathematical model was calibrated using experimental data from baboons administered a two-hour infusion of E coli and followed for a maximum of 28 days. Ensembles of parameters were generated using a Bayesian parallel tempering approach to produce model fits that could recreate experimental outcomes. Stepwise logistic regression identified seven model parameters as key determinants of mortality. Sensitivity analysis showed that parameters controlling the level of killer cell neutrophils affected the overall systemic damage of individuals. To evaluate rescue strategies and provide probabilistic predictions of their impact on mortality, time of onset, duration, and capture efficacy of an extracorporeal device that modulated neutrophil phenotype were explored. Our findings suggest that interventions aiming to modulate phenotypic composition are time sensitive. When introduced between 3–6 hours of infection for a 72 hour duration, the survivor population increased from 31% to 40–80%. Treatment efficacy quickly diminishes if not introduced within 15 hours of infection. Significant harm is possible with treatment durations ranging from 5–24 hours, which may reduce survival to 13%. In severe sepsis, an extracorporeal treatment which modulates CXCR-1/2 levels has therapeutic potential, but also potential for harm. Further development of the computational model will help guide optimal device development and determine which patient populations should be targeted by treatment